Singularity in self-energy and composite fermion excitations of interacting electrons 
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We propose that a composite fermion operator fia{2nia — 1) could liave coiierent excitations, where fi^ 
is the fermion operator for interacting electrons and rii^ is the number operator of the opposite spin. In the 
two-impurity Anderson model, it is found that the excitation of this composite fermion has a pseudogap in the 
Kondo regime, and has a finite spectral weight in the regime where the excitation of the regular fermion ficr has 
a pseudogap. In the latter regime, the self-energy of fia is found to be singular near Fermi energy. We argue that 
this composite fermion could develop a Fermi surface with Fermi liquid behaviors but "hidden" from charge 
excitations in lattice generalizations. We further illustrate that this type of excitations is essential in addressing 
the pseudogap state and unconventional superconductivity. 
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Introduction. The central issue in strongly correlated elec- 
tron systems is on how the strong Coulomb interaction leads 
to many fascinating phenomena such as magnetic quantum 
phase transitions (QPTs), Mott transitions, and unconven- 
tional superconductivity in strongly correlated electron sys- 
tems. In cuprate superconductors, a hotly debated issue is 
on the pseudogap state 111], as a doped Mott insulator, where 
the single-particle excitations are suppressed near the Fermi 
energy. Meanwhile, It is suggested the being in the vicinity 
of a Mott transition plays an important role in understand- 
ing the superconductivity in iron-based superconductors iQl- 
The "Kondo-breakdown" QPT in heavy fermions can also be 
viewed as a Mott transition in /-electrons Ql- The pseudo- 
gap physics, though not explicitly observed in the latter two 
systems, is expected to be present as well 101 ■ To understand 
the origin of a pseudogap therefore presents a key role in un- 
derstanding other emerging properties of correlated electrons, 
in particular, superconductivity. The major challenge, which 
prevents such a understanding, is the lack of proper theoret- 
ical methods to treat strong Coulomb interactions between d 
or / electrons in these materials, as captured by some funda- 
mental theoretical models such as the Hubbard model and the 
periodic Anderson lattice model. In this sense, the dynami- 
cal mean-field theory (DMFT) and its cluster extensions have 
made a great stride in this direction, by mapping the lattice 
problem into a self-consistent single (or a cluster) Anderson 
impurity model jsf]. In cluster-DMFT solutions to the Hub- 
bard model, a pseudogap is typically identified Further it 
has been shown that the pseudogap is associated with a sin- 
gularity in self-energies [7(1. However, there are no explicit 
explanations or detailed analyses on this behavior 

In this Letter, we provide such an analysis from a mini- 
mal two-impurity Anderson model iQ]. This model captures 
both the local quantum fluctuations such as the Kondo dy- 
namics, and the short-range spatial fluctuations, such as the 
inter-site spin or change correlations, which are also the key 
fluctuations in lattice systems. This model is subject to ex- 
act numerical solutions, such as the numerical renormaliza- 
tion group (NRG) method ^ as we adopt below. With the 
advantages of the NRG method in its low energy resolution 



and its dealing with real frequency directly, we obtain for the 
first time an accurate form for self-energies over the entire 
energy range for models with both local and nonlocal inter- 
actions. We show that the self-energy for Anderson orbitals 
becomes singular (with a pole) near the Fermi energy in cer- 
tain regimes controlled by inter-site spin exchange interaction 
(RKKY interaction). To address the origin of this singularity, 
we calculate the spectra of composite fermions, dia- = fianis 
and eic = fia{l - nm), where /^^ = + e^^ is the 
fermion operator for Anderson orbitals. The two composite 
fermion operators differentiate the single-particle excitations 
when the orbital is already single-occupied or vacant. We re- 
alize that when these two composite fermions have different 
behaviors, it is necessary to incorporate their anti -bonding op- 
erator fi„ = di„ - Cia = ftai'in^s - 1) to account for the 
single-particle excitations. We find that in the Kondo regime 
dominated by onsite Kondo coupling, with a Kondo resonance 
peak in fi^r, the excitations of fi^ have a pseudogap. In cer- 
tain regime dominated by RKKY interaction, where there is 
a pseudogap in excitations of fi^, the excitations of fic has 
finite spectral weight. We argue fi^, as a canonical fermion, 
could develop a Fermi surface with Fermi liquid behaviors but 
"hidden" from single charge excitations. The scattering from 
this composite fermion leads to a pseudogap (or singularity in 
self-energy) of fi^- We further show that a Cooper pair exci- 
tation is finite only when both excitations of fi^ and fi^ are fi- 
nite, implying that the mechanism of superconducting pairing 
may rely on the composite fermion excitations in interacting 
systems. 

Model and Method. We consider the Hamiltonian for the 
two-impurity Anderson model; 

kcr k(TZ 
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Here Cko- and fi„ (z = 1, 2) are electron operators conduction 
band and Anderson orbitals, with respective band dispersion 
Ek and energy level e/. The electrons on the local orbitals 
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experiences an onsite Coulomb repulsion U, while the non- 
interacting electron band hybridizes with local /-electrons on 
each site with the strength Vk- is the chemical potential. In 
applications to cluster-DMFT, fi„ are chosen from neighbor- 
ing sites in bipartite sublattices with Ck<j descrbing an effective 
fluctuating bath. In terms of the even or odd parity combina- 
tions of the local orbitals, /p=(e,o)<T = (./i<t ± /2(t)/a/2, the 
bath spectral spectrum Ae,o{uj) = Ek(^k /2A^L)|e*''''' ± 
g-tk r2 p^^-j^ — £!(; + *0"'"), is to be determined from a self- 
consistent procedure. In this case, the even and odd parity 
states correspond to the momentum points F = (0,0,..) and 
A/=(7r,7r,...) for a d-dimension lattice, respectively. 

We adopt the NRG method to solve this model. The cal- 
culation details can be found in Ref. [sl In particular, we 
carry out calculations of the self-energies for Anderson or- 
bitals I]pcr(a;) from the equation of motion, 

'^J^ai^) = w + /i-e/- A(w)-S;p^(w), 

SpaH = (([/p.,if/,,„d;/pU)/G/.p.(c.) . (2) 

As Hf^tnt = J2iU^fi't^fii' [fp<J,Hf_int]/U = dpcr = 

(/icr"-iCT ± f 2(7^20-) / V2 is a composite fermion operator. 
Fpa{uj) = {{[fpa,Hf,int]; f^„}} is usually calculated by 
a mean-field cutoff in analytical methods, but can be de- 
termined exactly in numerics, as firstly introduced for the 
single-impurity Anderson model | flol] . In NRG, the imagi- 
nary parts of dynamical quantities are directly calculated with 
the Lehmann representation. The real parts can be determined 
subsequently from the Kramers-Kronig relation. 

Results of self-energies. We start from a presumed form of 
the bath spectrum, Te.oioj) = — ImA(a;) = Tq for |a;| < D, 
where D is the bath electron bandwidth and is set as the en- 
ergy unit. This case has been studied earlier in Refs. [H] and 
d. We follow the same numerical procedure and adopt the 
same parameters Tq — 0.0457r, U = — 2e/ ~ 2 and /i = 
as in Ref. HI [as Case (i)]. As the generated RKKY interac- 
tion vanishes for this spectrum, we add an explicit intersite 
spin exchange term /Si • S2 to simulate the RKKY interac- 
tion effect. It is found that for I < Ic ~ 2.3T^, as the 
single-ion Kondo temperature, the ground state is a Kondo 
resonance state with finite spectrum Af{uj) = —lmGf{uj) at 
u! = but with reduced Fermi liquid temperature; for I > Ic, 
the ground state is dominated by the RKKY interaction with 
vanishing spectrum weight Af{uj) ^ iJ^ . At the quantum 
critical point (QCP) / = /c, various correlation functions are 
found to be divergent, including the staggered magnetic sus- 
ceptibility, the inter-site singlet Cooper pair correlation func- 
tion, and the current fluctuation between two sites III III12II . In 
Fig. 1, we show the results of self-energies for different val- 
ues of / [See also Fig. 3 for In the Kondo regime 
(/ < /c), it is found that the self-energy has an analytical 
form ReI]po-(i^) = U/2 — ujjZ [See inset of Fig. 1(a)] and 
ImSpo-(i^) ~ — at low energies, while Z — > approaching 
the QCR In the RKKY regime (/ > /c), where a pseudo- 
gap form for the single-particle spectra is identified, we find 
Rel]pcr(w) = Z' jijj with Z' ^ approaching the QCP, i.e.. 



the self energy has a pole at the Fermi energy. [We experi- 
ence 0/0 situation in determining ImI]po-('^) at — s- 0, whose 
exact form at low energies relies on numerical precision be- 
yond current calculation. From ReSpcrCi^) ^ we expect 
ImSpcr(i^) ~ const.]. Apparently, the development of a pseu- 
dogap is associated with this singularity in the self-energy, as 
both ReGpcr(<^) and YmGpaiuS) vanish at the energy where 
ReSpo-(i^) becomes infinite. 
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FIG. 1: (Color online) Real (a) and imaginary (b) parts of the self- 
energy for /-electrons for the case of V^^o(y^) = To with added 
/Si • S2 term. Here Tq = 0.0457r, U = -2ey = 2 and = 
0. Various quantities for different parity and spin channels are the 
same due to the symmetry. Additional particle-hole symmetry is also 
preserved. In the inset of panel (a), we show the zoom-in region near 
the Fermi energy in linear-linear scale to demonstrate the behavior of 
ReE in the Kondo regime. 

In Fig. 2, we further show the self-energies for another 
case with the bath spectrum determined from a realistic 3D 
tight-binding dispersion for conduction electrons. Here, it is 
found that re,o(w) have the formof ro(l=F'^) at low energies, 
and a finite antiferromagnetic RKKY interaction is perturba- 
tively generated [S]. Therefore, we turn off the explicit spin 
exchange term /Si • S2, and simply change the value of the 
hybridization constant Vk (with U fixed) to tune the relative 
strength between / and T^. Here the results are for [/ = 2, 
(j = —0.5, and /i = —0.2. Due to the lift of the symme- 
try between even and odd parities, there is no sharp transi- 
tion in this case. Instead, the system changes smoothly from 
a Kondo resonance state to an inter-impurity singlet state, 
respectively with finite and almost vanishing single-particle 
spectral weights near the Fermi energy. We also identify the 
singularity in the self-energy of /-electrons in certain range of 
the RKKY dominating regime. In this case, the singularity (as 
a pole) only exists in the even-parity channel near the Fermi 
energy. 

Composite fermion excitations. From Eq. (2), the emer- 
gence of a pole in the self-energy takes place when both 
ReG'pCT(a;) and ImG'pcr(a;) vanish while Fpcjiy) is finite. This 
motivates us to study the behaviors of separate composite 
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FIG. 2: (Color online) (a) Spectral function Af{u}), (b) real and 
(c) imaginary parts of the self-energy for different values of the hy- 
bridization constant V for a realistic band model. Here U = 2, 
ef — —0.5, and fi = —0.2. As V increases, the relative ratio be- 
tween the generated RKKY interaction / and decreases, and the 
system evolves from the RKKY regime to the Kondo regime. The 
solid and dashed lines represent the even and odd parity channels, 
respectively. The spin symmetry is preserved. The pole position in 
ReE for V = 0.15 is illustrated as example, with a (5-function-like 
peak in ImE. 



fermion operators dp^r = [dia ± d2<y]/V^ and Cpa = \e\a ± 
&ia\l\f^- [Notice that d^^ ^ fp^nps, similarly, for dko- in 
a lattice]. In the atomic limit (without the bath), / electrons 
only have the upper and lower Hubbard levels, which can be 
treated as the energy level for dia and Cia composite fermions, 
i.e., the /-electron is fractionaUzed into two flavors of com- 
posite fermions. As there is no hybridization between these 
two composite fermions, the single-particle excitation has a 
Mott gap, while there is also a pole in self-energy. 

In Fig. 3, we show the properties of the composite fermions 
as well as their hybridizations for the first case. In the 
Kondo regime (/ = 0.002), it is found that Im((dp„, e^^)) « 
Im((dp^,dt^)) w Im((/p^,/t^))/4attj = 0. Together with 
both vanishing Re(((ipcr, ej,^.)) and Kt{{dpa i d\,crl) ^ this im- 
plies that the hybridization between dp^r and Cp^r is the coher- 
ent part of the single-electron excitation /per, which behaves as 
a non-interaction electron with a renormalized energy level at 
the Fermi energy. Interestingly, their anti-bonding part fp^ = 
dpa — Spa has a pseudogap form, Aj{u)) = — ImGjr(w) ^ w^. 
It can be easily verified that /per is a canonical fermion, satis- 
fying the anti-commutation rule [/po-, /J/o-'l-f ~ ^pp'^au' [the 
same for the site basis]. In the RKKY regime (/ = 0.003), 
where Gf{uj) has a pseudogap, Af ^ w^, we find that Aj:{ijj) 
is finite near the Fermi energy. The latter resembles a "Kondo 
resonance" but with a smaller value for Ajr(O) [It is found to 
increase as V/U increases but its exact dependence is to be 
determined]. In this case, we find that lm{{dpcr, )) has an 



FIG. 3: (Color online) Imaginary and real parts of different com- 
posite fermion Green's functions. This is corresponding to the case 
in Fig. 1 for two values of J = 0.002 (a) and / = 0.003 (b). 
Lines with different colors represent different Green's functions: 
Gf = ((/p<T. /p<t)) (black), {{dpa,fla)) (red) [for comparison 
FpaUj/U (red dotted) with additional contributions from 7Si • S2 
term], {{dp^,dl^))^ (green), {{epa,el„)) (blue), {{dpa.el^)) (or- 
ange), and ((/pCT, Ha)) (cyan). -Im((dp<^, e|,„)) is negative in the 
energy range \u\ > C//2, and additionally near = for / = 0.003 
(not shown due to the log plot). Similarly, — Im{(dpo-, /po-)) is also 
negative in certain energy ranges. 
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FIG. 4: (Color online) The imaginary part of the Cooper pair cor- 
relation function Xac (green lines), the single-particle Green's func- 
tion Gf (black lines), and the Green's function Gj for the composite 
fermion fpa (red lines) for three representing values of RKKY inter- 
action in I < Ic, I ~ Ic, and / > Ic- 

opposite sign with lm{{dpa, d^pg.)), the opposite to the Kondo 
regime. In addition, Re((c?pcr, dpcr)) and Re((epcr, e^^)) take 
finite values of the order of This is a reminiscent of the 

atomic limit result, where, e.g., {{dpo- , dj,^)) = {npg)/{uj — 
£/ - U) but with lm{{dp„ , e^^)) = 0. Here Im((dpa, et^)) 
is rather finite at a; sa 0, serving as the difference between 
a pseudogap and a fully-opened Mott gap. The similar fea- 
tures can also be identified in the realistic band case, as shown 
in Fig. 2, but we find that away from single-occupancy, we 
need to redefine ( gp^ , 5pcr ) ~ ctdp„ ± Pcp^ with a//3 7^ 1 to 
capture the above duality (not shown). 

We proceed to study the superconducting instability. In 
Fig. 4, we show the Cooper pair correlation functions Xsc = 
((A,,,; AD), where Aj^ = fUfL+flfflv together with the 
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single particle spectra for f^^ and fj^^. We notice that the fi- 
nite pair excitations are closely related to the existence of both 
finite excitations of f^^ and /^^. In the particle-hole symmet- 
ric case, both A/ (0) and Aj(0) ai-e finite only at the QCP. But 
in the particle-hole asymmetric (or doped) case, they can be 
both finite in a range of tuning parameters. 

Discussions. As illustrated by the composite fermion exci- 
tations, the difference between a Kondo resonance state and 
a pseudogap state of fp„ can be understood as how the exci- 
tations of dprr and e.pa- are combined to form quasiparticles. 
In the Kondo regime, the elementary excitation is predomi- 
nantly in the form of the regular electron operator fpa- In the 
RKKY regime, it is rather in the form of an antibonding op- 



erator, fpo 



fpcr, as a canonical fermion, is ex- 



pected to develop a Fermi surface in the lattice generalization 
(or participate in forming the Fermi surface in heavy fermion 
problem). This also accounts for the Fermi liquid behaviors 
in the RKKY regime, with a finite Fermi liquid temperature 
vanishing at the QCP 11311 . But as the factor (2?!^^ — 1) « in 
the single-ocupancy limit, this composite Fermi surface is not 
sensitive to charge excitations involving single-particle exci- 
tations, such as current, electron tunneling. But it will con- 
tribute to spin dynamics and thermodynamics. In this sense, it 
can be dubbed as a "hidden"-Fermi surface. We notice that it 
is quite different from similar proposals from other theoretical 
aspects, such as the spinon Fermi surface. It does not rely on 
spin-charge separation and it is a Fermi liquid fixed point. We 
also notice that the composite fermion is in general itinerant, 
different from local moment pictures. 

On the other hand, the composite fermion operator /p^ 
can be treated as the original electron fermion operator com- 
bined with a bilinear electron operator capturing different 
particle-particle or particle-hole excitations, including charge 
and spin fluctuations. Such fluctuations have been treated as 
the bosonic fluctuations coupled to electrons. But such ap- 
proaches, emphasizing on low energy spin fluctuations, could 
never capture accurately the excitations of /p^- A bosonic 
fluctuation formalism could not produce a pole in the self- 
energy, but rather certain singular non-Fermi liquid form. In 
our calculations, all fluctuations, including Cooper pair fluc- 
tuations and current fluctuations are captured exactly and in- 
herentl y. W e notice a recent proposal of composite fermion 
theory |ll4f|, based on the Kotliar-Ruckenstein slave boson ap- 
proach, which can also capture the charge fluctuations and 
produce a pole in self-energy in terms of the hybridization 
between the composite fermion and the quasiparticles. Al- 
though we believe that there is a connection between their 
composite fermions with fpc here, a future study to estab- 
lish the exact relation is highly demanded. This composite 
fermion operator can be thought of as a dual representation 
of the original electron operator, i.e., the original Hamilto- 
nian can be represented by fpa only. The objective to intro- 
duce dp„ and Cpa is to eliminate the interaction term with 
quartic fermion operators. Therefore, as long as the exci- 
tations of dpa- and e^cr have different behaviors, it is neces- 
sary to incorporate fpa in an effective low-energy Hamilto- 



nian in terms of bilinear electron operators. We notice that, as 

fea = fea{nea + - 1) + flgfosfoa " feaflsfoa, VarioUS 

charge and spin fluctuations behave as a "hybridization" term 
llafecr- Howcvcr, duc to the nontrivial anti-commutation re- 
lation between fi^p and /^ct, it is not obvious how to construct 
such a hybridization model. 

We also find that the composite fermion does not have 
finite spectral weight in the local moment regime of the pseu- 
dogap Anderson impurity model, implying that the coherent 
excitation here is an inter-site interaction effect. How ex- 
actly the intersite spin exchange interaction leads to the sign 
change in lm{{dpa-, e^p„)), which directly differentiate the reg- 
ular fermion or the composite excitations, is to be studied in 
the future. 

In summary, we have proposed a canonical composite 
fermion, whose excitations are important to elucidate the sin- 
gularity in self-energy and the pseudogap state of the interact- 
ing electrons. 
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